home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / ddisna.z / ddisna
Encoding:
Text File  |  2002-10-03  |  4.1 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDDDDDIIIISSSSNNNNAAAA((((3333SSSS))))                                                          DDDDDDDDIIIISSSSNNNNAAAA((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DDISNA - compute the reciprocal condition numbers for the eigenvectors of
  10.      a real symmetric or complex Hermitian matrix or for the left or right
  11.      singular vectors of a general m-by-n matrix
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DDISNA( JOB, M, N, D, SEP, INFO )
  15.  
  16.          CHARACTER      JOB
  17.  
  18.          INTEGER        INFO, M, N
  19.  
  20.          DOUBLE         PRECISION D( * ), SEP( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      DDISNA computes the reciprocal condition numbers for the eigenvectors of
  37.      a real symmetric or complex Hermitian matrix or for the left or right
  38.      singular vectors of a general m-by-n matrix. The reciprocal condition
  39.      number is the 'gap' between the corresponding eigenvalue or singular
  40.      value and the nearest other one.
  41.  
  42.      The bound on the error, measured by angle in radians, in the I-th
  43.      computed vector is given by
  44.  
  45.             DLAMCH( 'E' ) * ( ANORM / SEP( I ) )
  46.  
  47.      where ANORM = 2-norm(A) = max( abs( D(j) ) ).  SEP(I) is not allowed to
  48.      be smaller than DLAMCH( 'E' )*ANORM in order to limit the size of the
  49.      error bound.
  50.  
  51.      DDISNA may also be used to compute error bounds for eigenvectors of the
  52.      generalized symmetric definite eigenproblem.
  53.  
  54.  
  55. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  56.      JOB     (input) CHARACTER*1
  57.              Specifies for which problem the reciprocal condition numbers
  58.              should be computed:
  59.              = 'E':  the eigenvectors of a symmetric/Hermitian matrix;
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDDDDDIIIISSSSNNNNAAAA((((3333SSSS))))                                                          DDDDDDDDIIIISSSSNNNNAAAA((((3333SSSS))))
  71.  
  72.  
  73.  
  74.              = 'L':  the left singular vectors of a general matrix;
  75.              = 'R':  the right singular vectors of a general matrix.
  76.  
  77.      M       (input) INTEGER
  78.              The number of rows of the matrix. M >= 0.
  79.  
  80.      N       (input) INTEGER
  81.              If JOB = 'L' or 'R', the number of columns of the matrix, in
  82.              which case N >= 0. Ignored if JOB = 'E'.
  83.  
  84.      D       (input) DOUBLE PRECISION array, dimension (M) if JOB = 'E'
  85.              dimension (min(M,N)) if JOB = 'L' or 'R' The eigenvalues (if JOB
  86.              = 'E') or singular values (if JOB = order. If singular values,
  87.              they must be non-negative.
  88.  
  89.      SEP     (output) DOUBLE PRECISION array, dimension (M) if JOB = 'E'
  90.              dimension (min(M,N)) if JOB = 'L' or 'R' The reciprocal condition
  91.              numbers of the vectors.
  92.  
  93.      INFO    (output) INTEGER
  94.              = 0:  successful exit.
  95.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  96.  
  97. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  98.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  99.  
  100.      This man page is available only online.
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.